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Abstract 

Belle report data on T(55) — >■ T(lS,2S,3S)ir + -K~ and x&(l-P, 2P)ir + -K~ ; they observe 
peaks in T-k^ and consistent with J p = 1 + . They interpret the peaks as molecular 

states Z;,(10608) and Z^(10653). Their masses are just above BB* and B*B* thresholds at 
10604.6 and 10650.2 MeV. An explanation in terms of cusps at these thresholds is presented 
' here. The product of the rising phase space for BB* and B*B* with the cusps creates peaks 

a few MeV higher in BB* and B*B*, and these peaks can de-excite to 7r + 7r~T(15, 25, 3S) 
and 7r + Tr-Xb(lP,2P). 
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Belle P report peaks in 7r ± Y(lS', 25*, 3S) and 7r Xf,(lP, 2P) with masses and widths: 



T^i M[Z 6 (10610)] = 10608.4 ±2.0 MeV T = 15.6 ± 2.5 MeV (1) 

M[Z fe (10650)] = 10653.2 ± 1.5 MeV T = 14.4 ± 3.2MeV. (2) 

Angular distributions are consistent with J p = 1 + in both cases. 

Thresholds for B*B and B*B* are at 10604.6 MeV and 10650.2 MeV and have the same mass 
, separation as the two peaks within experimental errors. This suggests strongly that the 
peaks are associated with those thresholds. There is an explanation how this can arise via cusps 
at the B*B and B*B* thresholds. 



Consider B*B* elastic scattering first. The imaginary part of the T-matrix for B*B* —s- B*B* 
S-wave scattering is 

ImT oc g 2 p(s), (3) 

where p is the phase space for B*B* and g is the coupling constant to this channel. Analyticity 
requires a real part to the amplitude [2] 

S: RerK Ip/™^^£>M, (4) 



7T Js thr S — S 



where s t hr is mass squared at threshold; P denotes the principal value integral. For B*B*, the 
S-wave phase space is 

2k 

p(s') = -=FF{J), (5) 



where k is the momentum of each B* in the B*B* rest frame and FF(s ; ) is a form factor. There 
are further form factors for T(55) — > B*B*, but these change little over the cusps and will be 
neglected. It is important to be aware that this cusp effect is not an 'optional extra'. It is a 
basic requirement at the opening of any new channel, but has been ignored in the majority of 
calculations of 'molecules'. It is also known as the 'coupled channel effect'. 
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Figure 1: One diagram for T(5S) ->■ n + B*B* followed by de-exitation B*B* ->■ 7r T(nS). Full 
lines represent 6 and 6, dashed lines light quarks. 



Let us now turn to the process T(5S) — > 7i + n^T(nS) shown in Fig. 1. It goes via re- 
arrangement of the four quarks in B*B*. This process has the same analytic structure as 
B*B* — > B*B*, since that amplitude obviously factorises. Let us consider first B*B* — > T(nS)n 
via 7r exchange. Heavier exchanges will be discussed later. The momentum transfer in the pro- 
duction is quite large and varies little over the cusp, though the phase space for the intermediate 
B* is included. For free B*, pion decay is forbidden by the small mass difference between B* 
and B. However, for the triangle graph shown in Fig. 1, intermediate B* are off-shell, allowing 
it exchange in B*B scattering at the right-hand vertex. Suppose the form factor is a Gaussian 

FF = exp (-k 2 R 2 /3) (6) 

with R — 1.41 fm, corresponding to the mass of the pion, though the exponential may be an 
approximation. 
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Figure 2: ImT, ReT or T/k (dotted) and |T| 2 (full curve) for (a) B*B ->• B*B and (b) B*B* ->■ 
B*B*; Int stands for Intensity. On Fig. 2(a), ReT is also shown by a full curve. 



Fig. 2(b) shows ImT as the dashed curve, normalised to 1 at its peak. Then ReT , shown as 
the dotted curve with the same normalisation, has a sharp cusp exactly at the B*B* threshold. 
ImT peaks at a mass of 10656 MeV. ReT goes negative just above this mass. For higher masses, 
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there is repulsion between B* and B*; this will keep them apart and suppress de-excitation to 
TiT(nS), where n = 1—3. As a result, the intensity \T\ 2 for T(5S f ) —> n + n^T(nS) is proportional 
to the full curve. It rises at the threshold 10650.2 MeV to the scattering length given by the 
cusp and switches off at 10658.7 MeV. Folding in the 5.2 MeV experimental resolution quoted by 
Belle gives a somewhat rounded line-shape with a full-width at half-height of 10.0 MeV, rather 
less than the Belle value of 14.4 ± 3.2 MeV. 

Consider next B*B scattering via a diagram like Fig. 1 except for a B~ running from left- 
to right-hand vertices. This requires a P-wave interaction and 



assuming a centrifugal barrier with the same radius as the exponential; this is an approximation 
where the centrifugal barrier is replaced by an equivalent square barrier. It is necessary to 
preserve the P-wave form for the amplitude at threshold and the dispersion integral becomes 



Fig. 2(a) shows the corresponding line-shape for KeT/k as the dotted curve; as on Fig. 2(b), 
it peaks at threshold. When it is multiplied by the factor k in Eq. (8), the line-shapes for ReT 
and \T\ 2 are shown by the full curves and ImT is shown by the dashed curve. In this case, 
the threshold is at 10604.6 MeV and the B*B interaction becomes repulsive at 10613.9 MeV. 
The peak is centred at 10609.2 MeV and its full width becomes 9.5 MeV when folded with 
experimental resolution, compared with the experimental value 15.6 ± 2.5 MeV. 

There are two key points in favour of the cusp mechanism. The first is that it predicts 
almost the same widths for Z&(10610) and -2^(10650), in agreement with data. In the molecular 
picture, these widths have no obvious relation. Secondly, the measured branching ratios of 
T(1S, 2S, 3S)7in sum to ~ 2% of the total width of T(10860), hence accounting for only a 2 
MeV width for each peak [3]. There is a factor 7 discrepancy with experimental values. 

Let us now consider heavier exchanges, for example due to a. A full calculation is plagued 
by unknowns. Firstly, it is unclear whether the a will couple via the a pole (as in BES II data 
for J/ty — > U7T7T [4j) or as in elastic scattering, where the Adler zero appears in the numerator 
of the amplitude. Next, it is necessary to fold the radius of interaction (as a function of a mass) 
with the folded radius of two B* (or BB*). A calculation gives an RMS radius of ~ 0.65 fm 
after this folding procedure. The cut-off at high mass varies with the folded mass and is 
rounded off both below and above the sharp cut-off as shown below in Fig. 3. However, the key 
point is that, classically, the cross section for each de-excitation process, e.g. B*B* — > T(nS)ii, 
is proportional to the mean < R 2 >. As a result 7r exchange plays a strong role. 

This procedure will be reversed here, deducing the RMS radius from the widths of ^(10610) 
and Zfe(10650) and asking whether this is reasonable. Even then, there is a minor problem 
that Belle have assumed these two peaks are resonant and included interferences between them 
and have also included interference with a non-resonant background. The two most prominent 
signals in their data are for Z;,(10610) and Zb(10650) — > T(2S, 3S)7m, see their Fig. 5. Reading 
data directly from these graphs, widths come out 17% higher than their values. This is a small 
discrepancy, but systematic and probably arises from their inclusion of interferences between 
the two peaks. Using my widths, V < R 2 > comes out to be 0.95 fm. In view of the factor R 2 
favouring ir exchange, this appears reasonable within ± 0.1 fm, i.e. ~ 10%. 



p(B*B) = k 3 



exp(-k 2 R 2 )/(l + k 2 R 2 ) 



(7) 
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Figure 3: Line-shapes of the cusps using a radius of 0.95 fm, adjusted to give an optimum fit to 
Belle data. 



Fig. 3 shows resulting line-shapes for this radius of interaction. The sharp cut-offs of Fig. 
2 are slightly rounded off, but are still distinctive. With further data of higher statistics (and 
perhaps improved mass resolution), this line-shape is open to experimental test. In work con- 
current with the present paper, Danilkin, Orlovsky and Simonov discuss the general formalism 
of heavy quarkonia [5] and show results for Z&(10610) and -Zfe(10650); they use a smaller radius 
of interaction and therefore find wider cusps which are inconsistent with the data. 

It is possible that the Zf, states are actually resonant. Cleven, Guo, Hanhart and Mei/3ner 
assume the molecular state is resonant and close to threshold; they adjust coupling parameters 
to reproduce the data on hb(2P) final states [6]. They include the cusp effect and arrive at quite 
different line-shapes to those presented here. Their essential feature is a sharp cusp right at 
threshold. 
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Figure 4: Argand diagrams for (a) B*B — > B*B, (b) B*B* — > B*B*; numbers are masses in 
GeV. 



Figs. 4(a) and (b) show Argand diagrams for both Zj, derived from Fig. 3; masses are marked 
in GeV. In Fig. 4(b) there is no pole, because of the s-dependence around the curve. Above 
a mass where the amplitude goes negative, the variation with s slows down rapidly. However, 
there is a possibility that meson exchanges will dynamically generate a resonant pole. The 
evaluation of the cusp made here is exactly equivalent to evaluating the first order loop diagram 
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for meson exchanges. Solving the Bethe-Salpeter equation would include loop diagrams of all 
orders and could dynamically generate a resonance. The more favourable situation is for the 
loop driven by BB*, because of the P-wave factor at threshold. For Fig. 4(a), if the data are 
fitted with a Breit-Wigner resonance with s-dependent width, there is a pole at 10.623 ± 0.005 
MeV. It has a full width at half-height of 44 MeV and is centred at 10632 MeV. Note that the 
loop is really a property of the BB* channel. If there is an isospin 1 = Xbi(P) state near this 
threshold, it is likely to be attracted towards the B* threshold. There is also the possibility that 
the Zb(10610) becomes an exotic resonance, but that needs to be proved. 

There is one detail concerning decays to 7i + 7i~T(nS). For this case, both 7r are produced in 
a relative S-wave and are fitted with the a amplitude. To form the final states 7r + 7r~xt,(lP, 2P) 
observed by Belle, one pion must be in a P-wave. This is consistent with the observation of 
decays to <tl=i, where L is the orbital angular momentum of the a with respect to Xb states. 
A further detail is that pion production in the step B*B* or B*B — > -n\b is likely to flip the 
spin of the Xb final state. This can account for a phase difference with respect to 7i + 7i^T(nS) of 
~ 180°, as recorded by Belle. Because of nodes in radial wave functions, significant differences 
are to be expected between Dalitz plots for 7r + 7rT(nS') with different n values. There has been 
earlier discussion of a possible hybrid in the mass range 10.4-10.8 GeV on the basis of variations 
of Dalitz plot for Ttttt with different n values [7] . 

In conclusion, threshold cusps at B*B and B*B* thresholds are capable of explaining the 
peaks observed by Belle. The predicted line-shapes are distinctive, and with higher statistics 
can be used to test the present hypothesis. The width predicted for these cusps requires some 
exchange of heavier mesons than pions between B* and both B and B*. The cusp mechanism 
has the virtue of predicting closely similar widths for both Zj, peaks. The experimental width 
of the peaks is a factor 7 larger than that for Z b — > T(1S, 2S, 3S)7m, but is reproduced by the 
cusp mechanism with a reasonable radius of interaction. 
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